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Abstract: Undulatory locomotion, as seen in the nematode Caenorhabditis ele- 
gans, is a common swimming gait of organisms in the low Reynolds number regime, 
where viscous forces are dominant. While the nematode's motility is expected to be 
a strong function of its material properties, measurements remain scarce. Here, the 
swimming behavior of C. elegans are investigated in experiments and in a simple 
model. Experiments reveal that nematodes swim in a periodic fashion and generate 
traveling waves which decay from head to tail. The model is able to capture the 
experiments' main features and is used to estimate the nematode's Young's modu- 
lus E and tissue viscosity rj. For wild-type C. elegans, we find E ~ 3.77 kPa and 
n ~ —860 Pa-s; values of n for live C. elegans are negative because the tissue is 
generating rather than dissipating energy. Results show that material properties are 
sensitive to changes in muscle functional properties, and are useful quantitative tools 
with which to more accurately describe new and existing muscle mutants. 



I. INTRODUCTION 

Motility analysis of model organisms, such as the nematode Caenorhabditis elegans (C. 
elegans), is of great scientific and practical interest. It can provide, for example, a powerful 
tool for the analysis of genetic diseases in humans such as muscular dystrophy (MD) [H 121 [3] 
since C. elegans have muscle cells that are highly similar in both anatomy and molecular 
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makeup to vertebrate skeletal muscles [HE]. Due to the nematode's small size (L « 1 mm), 
the motility of C. elegans swimming in a simple, Newtonian fluid is usually investigated in 
the low Reynolds numbers (-Re) regime, where linear viscous forces dominate over nonlinear 
inertial forces (6j Ej. At low Re, locomotion results from non-reciprocal deformations to 
break time- reversal symmetry [8j; this is the so-called "scallop theorem" [9]. Experimental 
observations have shown that motility of swimming nematodes including C. elegans results 
from the propagation of bending waves along the nematode's body length [101 [HI [12]. These 
waves consist of alternating phases of dorsal and ventral muscle contractions driven by the 
neuromuscular activity of muscle cells. While it is generally accepted that during locomotion 
the nematode's tissues obey a viscoelastic reaction [131 D3J US], quantitative data on C. 
elegans' material properties such as tissue viscosity and Young's modulus remain largely 
unexplored. 

Motility behavior of C. elegans is a strong function of its body material properties. 
Recent investigations have provided valuable data on C. elegans* motility, such as velocity, 
bending frequency, and body wavelength [TTJ [TJl EH EH EH EE]- However, only recently 
have the nematode's material properties been probed using piezoresistive cantilevers [19]. 
Such invasive measurements provided Young's modulus values of the C. elegans ' cuticle on 
the order of 400 MPa; this value is closer to stiff rubber than to soft tissues. 

In this paper, we investigate the motility of C. elegans in both experiments and in a 
model in order to estimate the nematode's material properties. Experiments show that 
nematodes swim in a highly periodic fashion and generate traveling waves which decay from 
head to tail. A dynamic model is proposed based on force and moment (torque) balance. A 
simplified version of the model is able to capture the main features of the experiments such 
as the traveling waves and their decay. The model is used to estimate both the Young's 
modulus and tissue viscosity of C. elegans. Such estimates are used to characterize motility 
phenotypes of healthy nematodes and mutants carrying muscular dystrophy (MD). 

II. EXPERIMENTAL METHODS 

Experiments are performed by imaging C. elegans using standard microscopy and a high- 
speed camera at 125 frames per second. We focus our analysis on forward swimming in 
shallow channels to minimize three-dimensional motion. Channels are machined in acrylic 
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and are 1.5 mm wide and 500 jum deep; they are sealed with a thin (0.13 mm) cover glass. 
Channels are filled with an aqueous solution of M9 buffer [20], which contains 5 to 10 
nematodes. The buffer viscosity p and density p are 1.1 cP and 1.0 g/cm 3 , respectively. 
Under such conditions, the Reynolds number, defined as Re = pUL/p, is less than unity, 
where U and L are the nematode's swimming speed and length, respectively. 

In Fig. [TJa), we display nematode tracking data over multiple bending cycles for a healthy, 
wild-type nematode. Results show that the nematode swims with an average speed < U >= 
0.45 mm/s and with a beating pattern of period T = 0.46 s. This periodic behavior is 
also qualitatively observed in the motion of the nematode tail (Fig. [T£i); see also Video 1 
(Supplementary Material). Under such conditions, Re ~ 0.4. Snapshots of the nematode 
skeletons over one beating cycle (Fig. [T|d) reveal an envelope of well-confined body postures 
with a wavelength of approximately 1 mm, which corresponds nearly to the nematode's 
body length. The displacement amplitudes at the head and tail are similar with 465 pm 
and 400 /im, respectively. However, the amplitudes of the curvature at the head and tail 
differ sharply with approximately 6.07 mm -1 and 2.21 mm -1 , respectively. The tail/head 
curvature ratio of 0.36 suggests that the bending motion is initiated at the head. 

Extensive genetic analysis in C. elegans has identified numerous mutations affecting ne- 
matode motility. One such mutant, dys-1, encodes a homolog of the human dystrophin 
protein, which is mutated in Duchenne's and Becker's muscular dystrophy (MD). Using qual- 
itative observation, dys-1 mutants have an extremely subtle movement defect [21] . which 
includes slightly exaggerated head bending and time-dependent decay in movement. The 
quantitative imaging platform presented here is able to robustly differentiate between wild 
type and dys-1 mutants, as shown in Fig. [ljc) and Video 2 (Supplementary Material). 

Results show that the dys-1 mutant swims with an average speed < U >= 0.17 mm/s and 
Re ~ 0.15; both values are significantly smaller than for the wild-type nematode. Although 
the dys-1 mutant suffers from severe motility defects [21], it still moves in a highly periodic 
fashion with T = 0.63 s. Snapshots of nematode skeletons over one beating cycle (Fig. [Tji) 
also reveal an envelope of well-confined body postures with a wavelength corresponding to 
the nematode's body length. The dys-1 mutant exhibits a tail/head curvature ratio of 0.23, 
which is similar to the value found for wild type nematodes. However, the corresponding 
displacement amplitudes of the mutant are much smaller than those observed for the wild- 
type (Fig. fib). The displacement amplitudes at the head and tail are 330 pm and 155 pm, 
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respectively. This observation suggests that the dys-1 mutant body, and in particular the 
tail, are becoming inactive as the bending motion is not able to deliver as much body 
displacement. 

To further characterize the motility of C. elegans, we measure the curvature n(s,t) = 
d(f>/ds along the nematode's body (Fig. [2^,). Here, <p is the angle made by the tangent to 
the x-axis at each point along the centerline and s is the arc-length coordinate spanning 
the nematode's head (s = 0) to its tail (s = L). The spatio-temporal evolution of k for a 
swimming nematode is shown in Fig. [2^a). Approximately 6 bending cycles are illustrated 
and curvature values are color-coded; red and blue represent positive and negative values of 
k, respectively. The y-axis in Fig. [2^a) corresponds to the non-dimensional body position 
s/L. The contour plot shows the existence of highly periodic, well-defined diagonally oriented 
lines. These diagonal lines are characteristic of bending waves, which propagate in time along 
the body length. Note that as the wave travels along the nematode body, the magnitude of k, 
decays from head to tail. Such behavior contrasts sharply with that observed for undulatory 
swimmers of the inertial regime (e.g. eel, lamprey), where amplitudes of body displacement 
grow instead from head to tail [22j E3]. 

The body bending frequency (/) is obtained from the one-dimensional Fast Fourier Trans- 
form (FFT) of the curvature field k at multiple body positions s/L (Fig. ^f>). Here, the 
body bending frequency is defined as / = oj/lit. The angular frequency uj is calculated by 
first extracting multiple lines from the curvature field at distinct body positions s/L, and 
then computing the one-dimensional FFT. The wave speed c is extracted from the slope of 
the curvature k propagating along the nematode's body; the wavelength A is computed from 
the expression A = c/f. A single frequency peak / = 2.17 ± 0.18 Hz (n = 25) is found in 
the Fourier spectrum, where n is the number of nematodes. This single peak is irrespective 
of body position and corresponds to a wave speed c = 2.14 ± 0.16 mm/s. 

III. MATHEMATICAL METHODS 

The swimming motion of C. elegans is modeled as a slender body in the limit of low 
Reynolds numbers (Re) [2U 1251 I2S]- This model is later used to estimate the material 
properties of C. elegans. We assume that the nematode is inextensible [26]; the uncertainty 
in the measured body lengths is less than 3%. The nematode's motion is restricted to the 
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xy-plane and is described in terms of its center-line y(s, t), where s is the arc-length along the 
filament [27]. The swimming C. elegans experiences no net total force or torque (moments) 
such that, in the limit of low Re, the dynamic equations of motion are 

OF 

— = C t ut + C n u n} (1) 
dM 

— = -[F y cos(4>) -F,sin(</>)]. (2) 

In Eq. F(s,t) is the internal force in the nematode, Cj is the drag coefficient experi- 
enced by the nematode, ui is the nematode velocity, and the subscripts t and n correspond 
to the tangent and normal directions, respectively. The drag coefficients Ct and C n are ob- 
tained from slender body theory [2E]. Due to the finite confinement of nematodes between 
parallel walls, corrections for wall effects on the resistive coefficients are estimated for slender 
cylinders [3 EH]. 

In Eq. g, M = M p + M a , where M p is a passive moment and M a is an active moment 
generated by the muscles of the nematode; the active and passive moments are parts of a 
total internal moment [T3J [H] . The passive moment is given by the Voigt model [15] such 
that M p = EIk + i] p I(dn/dt), where / is the second moment of inertia of the nematode 
cross section. The Voigt model is one of the simplest models for muscle and is extensively 
used in the literature [29]. Qualitatively, the elastic part of the Voigt model is represented 
by a spring of stiffness E while the dissipative part of the Voigt model is represented by 
a dashpot filled with a fluid of viscosity r] p (Fig. [3]). Here, we assume two homogeneous 
effective material properties, namely (i) a constant Young's modulus E and (ii) a constant 
tissue viscosity rj p . 

The active moment generated by the muscle is given by M a = —(ElK a + i] a IdK,/dt), 
where K a is a space and time dependent preferred curvature produced by the muscles of 
the nematode and r] a is a positive constant [30]. A simple form for n a can be obtained by 
assuming that n a is a sinusoidal function of time with an amplitude that decreases from 
the nematode's head to its tail (see Appendix). Note that if r\ = r\ p — r\ a > 0, there is net 
dissipation of energy in the tissue; conversely, if 77 = r] p — r] a < 0, there is net generation 
of energy in the tissue. Experiments have shown that the force generated by active muscle 
decreases with increasing velocity of shortening [29J, so that the force- velocity curve for 
active muscle has a negative slope (Fig. [3]). Such negative viscosity has been derived by a 
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mathematical analysis of the kinetics of the mechano-chemical reactions in the cross-bridge 
cycle of active muscles |30j . For live nematodes, we expect r\ = r\ v — r\ a < because the net 
energy produced in the (muscle) tissue is needed to overcome the drag from the surrounding 
fluid. 

Equations ([!]) and ^ are simplified by noting that the nematode moves primarily in the 
x-direction (see Video 1, S.M.) and that the deflections of its centerline from the x-axis are 
small. In such case, s ~ x and cos(</>) ~ 1. This results in a linearized set of equations given 
by 

Differentiating Eq. Q with respect to x and combining with Eq. ([3]), we obtain 

d 2 M dy 

-dx^ + Cn d-t = °- (5) 

Substituting for M(x,t) in terms of n(x,t) and its time derivative yields a bi-harmonic 
equation for the displacement y(x, t) of the type 

d A y dy 

which can be solved analytically for appropriate boundary conditions, where £ is a con- 
stant that depends on the nematode's material properties and the fluid drag coefficient (see 
Appendix). 

The boundary conditions are such that both the force and moment at the nematode's 
head and tail are equal to zero. That is, F y (0,t) = F y (L,t) = and M(0, t) = M(L,t) = 0. 
Note that the zero moment boundary conditions at the head and tail imply that EIk(0, t) + 
r]I(dK(0,t)/dt) = ElK a (0,t) and EIk(L, t) + -qI(dK(L,t)/dt) = EIn a (L, t). 

Experiments show that the curvature k has non-zero amplitudes (Fig. [2^l) both at the head 
(x = 0) and the tail (x = L). In order to capture this observation, we assume that K a (x, t) is a 
sinusoidal wave with decreasing amplitude of the form n a (x, t) = Q cosut + QiX cos(u;t — B), 
where Qq, Q\ and B are inferred from the experiments. Note that if the curvature amplitude 
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at the head is larger than that at the tail, then the nematode swims forward. Conversely, 
if the curvature amplitude is smaller at the head than at the tail, the nematode swims 
backward; if the amplitudes are equal at the head and tail then it remains stationary. We 
note, however, that other forms of the preferred curvature n a are possible and could replicate 
the behavior seen in experiments. 

IV. RESULTS & DISCUSSION 

Equation (|6| is solved for the displacement y(x, t) in order to obtain the curvature 
K(x,t) = d 2 y/dx 2 . The solution for y(x,t) is a superposition of four traveling waves of 
the general form A;exp(— /3xcosPj) cos(/3a;sinPj — cot — (pi) where (3 = (C n tu / 'K^) 1 ^ and 
Pi is a function of the phase angle tp. The amplitude A4 and phase fa are constants to 
be determined by enforcing the boundary conditions discussed above (see Appendix). The 
solution reveals both the traveling bending wave and the characteristic decay in k, as seen in 
experiments. Note that our formulation does not assume a wave-functional form for n(x,t). 
Rather the wave is obtained as part of the solution. 

Next, the curvature amplitude \k(x)\ predicted by the model is fitted to those obtained 
from experiment to estimate the bending modulus Kf, = I\J E 2 + uj 2 rf and the phase angle 
ip = ta,n~ 1 (r]Lj)/E. The nematode is assumed to be a hollow, cylindrical shell [T^l EI] such 
that I = ?r((r m + t/2) 4 - (r m -t/2) 4 )/4, where the mean nematode radius r m « 35 /im and 
the cuticle thickness t ~ 0.5 /im [32]. For the population of wild-type C. elegans tested here 
(n = 25), the best fit values are K b = 4.19 x i(T 16 ±0.49 x KT 16 Nm 2 and ip = -45.3° ±3.0°. 
In Fig. [4j the experimental values of \k(x)\ along the body of a wild-type C. elegans are 
displayed together with theoretical values of |«(a;)|, which are obtained by using the best 
fit value of the bending modulus for this nematode and by changing the phase angle from 
ip = 0° to ip = —90°. Figure [1] shows that the model is able to capture the decay in \k\ as a 
function of body length and the nematode's viscoelastic behavior. 

The values of Young's modulus E and tissue viscosity 77 can now be estimated based on 
the values of ip and K b discussed above. Results show that, for the wild-type nematodes, 
E = 3.77 ± 0.62 kPa and r] = -860.2 ± 99.4 Pa-s. The estimated value of E lies in the 
range of values of tissue elasticity measured for isolated brain (0.1 — 1 kPa) and muscle 
cells (8 — 17 kPa) [33]. The values of rj for live C. elegans are negative because the tissue 
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is generating rather than dissipating energy [301 EH ESI EE]. We note, however, that the 
absolute values of tissue viscosity \rj\ are within the range (10 2 — 10 4 Pa-s) measured for 
living cells [371 EH]- In order to determine whether the nematode's material properties can 
be extracted reliably from shape measurements alone, experiments in solutions of different 
viscosities were conducted [39]. The inferred Young's modulus and effective tissue viscosity 
remain constant for up to a 5-fold increase in the surrounding fluid viscosity (or mechanical 
load) . 

The nematode's curvature n(x, t) is now determined from the estimated values of E and 
r\. Figure [2|c) shows a typical curvature k(x, t) contour plot obtained from the solution of 
the above equations using the estimated values of E and rj. While the influence of non- 
linearities is neglected for the scope of the present paper, the analytical results show that 
our linearized model, while not perfect, is nevertheless able to capture the main features 
observed in experiments (Fig. |2]i). 

Next, the method described above is used to quantify motility phenotypes of three distinct 
mutant MD strains (see Table 1 in Supplementary Material): one with a well-characterized 
muscle defect (dys-1 ;hlh-l); one with a qualitatively subtle movement defect (dys-1), and 
one mutant that has never been characterized with regards to motility phenotypes but 
is homologous to a human gene that causes a form of muscular dystrophy expressed in 
nematode muscle (fer-1). Note that while both fer-1 and dys-1 genes are expressed in C. 
elegans muscle, they exhibit little, if any, change in whole nematode motility under standard 
lab assays [ID] . 

Figure [5] displays results of both kinematics (|5^) and tissue material properties (]5b) for all 
nematodes investigated here. Quantitative results are summarized in Table 1 (Supplemen- 
tary Material). We find that all three mutants exhibit significant changes in both motility 
kinematics and tissue properties. For example, fer-1 mutants exhibit defects in motility 
kinematics which are not found with standard assays [3D]. Specifically, both the maximum 
amount of body curvature attained in fer-1 mutants is increased by ~ 5%, and the rate 
of curvature decay along the body is increased by ~ 5% (Fig. [5^,). This data show that 
fer-1 (hc24) mutants exhibit small yet noteworthy defects in whole nematode motility and 
exhibit an uncoordinated (unc) phenotype. In comparison, kinematics data on dys-l;hlh- 
1 show that body curvature at the head of such mutant nematodes increases by ~ 70% 
compared to wild-type nematodes, while the rate of decay along the body is increased by 
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approximately ~ 40%. These results are useful to quantify the paralysis seen earlier in the 
tail motion of such MD mutants (Fig. [Tj^) . 

The Young's modulus (E) and the absolute values of tissue viscosity (\t)\) of wild-type 
and mutant strains are shown in Fig. [5](b). Results show that mutants have lower values 
of E when compared to wild-type nematodes. In other words, dys-1, dys-l;hlh-l, and fer-1 
mutants C. elegans are softer than their wild- type counterpart. The values of \r)\ of fer-1 
mutants are similar to wild-type nematodes, within experimental error. However, the values 
of | T] | for dys-1 mutants are lower than wild- type C. elegans. Since muscle fibers are known 
to exhibit visible damage for dys-1 ;hlh-l mutants [21 j, we hypothesize that the deterioration 
of muscle fibers may be responsible for the lower values of E and \rj\ found for dys-1 and 
dys-1 ;hlh-l mutants. 



V. CONCLUSION 



In summary, we characterize the swimming behavior of C. elegans at low Re. Results 
show a distinct periodic swimming behavior with a traveling wave that decays from the 
nematode's head to tail. By coupling experiments with a linearized model based on force 
and torque balance, we are able to estimate, non-invasively, the nematode's tissue material 
properties such as Young's modulus (E) and viscosity (77) as well as bending modulus (K b ). 
Results show that C. elegans behaves effectively as a viscoelastic material with E ~ 3.77 kPa, 
\r)\ ~ 860.2 Pa-s, and Kb ~ 4.19 x 10~ 16 Nm 2 . In particular, the estimated values of E are 
much closer to biological tissues than previously reported values obtained using piezoresistive 
cantilevers [19]. We demonstrate that the methods presented here may be used, for example, 
to quantify motility phenotypes and tissue properties associated with muscular dystrophy 
mutations in C. elegans. Overall, by combining kinematic data with a linearized model, 
we are able to provide a robust and highly quantitative phenotyping tool for analysis of 
C. elegans motility, kinematics, and tissue mechanical properties. Given the rapid non- 
invasive optical nature of this method, it may provide an ideal platform for genetic and 
small molecule screening applications aimed at correcting phenotypes of mutant nematodes. 
Our method also sheds new light on our understanding of muscle function, physiology, and 
animal locomotion in general. 
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APPENDIX: 

In this appendix, we detail the model for the motion of the nematode C. elegans. The 
nematode is modeled as a slender filament at low Reynolds numbers [22] • In our experiments, 
the uncertainty in the measured body lengths is less than 3% and we assume inextensibility. 
The nematode's motion is described in terms of its center-line y(s,t), where s is the arc- 
length along the filament and t is time [27]. We assume that the nematode moves in the 
xy-plane. The swimming C. elegans experiences no net total force or torque (moments) such 
that, in the limit of low Re, the equations of motion are 

OF 

— = C t U t + C n U n , (A.l) 

°^e z = -t x F, (A. 2) 

where t(s,t) = dy/ds is the tangent vector to the center-line, F(s,t) is the internal force, 
and M(s,t)e z = M(s,t) is the internal moment consisting of a passive and active part 
[T3| HI) . Tangential and normal velocities are respectively given by vit = (dy/dt ■ t)t and 
u n = (l—i®t)dy/dt. The drag coefficients, C t and C n , are obtained from slender body theory 
[2T?] . Due to the finite confinement of nematodes between the parallel walls, corrections for 
wall effects on the resistive coefficients are estimated for slender cylinders [7]. The local 
body position y, the velocity at any body position dy/dt, and the tangent vector t are all 
experimentally measured. 

The constitutive relation for the moment M(s, t) in our inextensible filament is assumed 
to be given by 

M = M p + M a , (A.3) 

where M p (s,t) is a passive moment and M a (s,t) is an active moment generated by the 
muscles of the nematode. The passive moment is given by a viscoelastic Voigt model [15] 



M p = ElK + ri P I-^, (A.4) 



where k(s, t) is the curvature along the nematode. Here, we assume two homogeneous 
effective material properties, namely (i) a constant Young's modulus E and (ii) a constant 
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tissue viscosity rj p . The active moment generated by the muscle is assumed to be given by: 

M a = -(ElK a + Va I—), (A.5) 

where n a = Qocosut + Qiscos(ut — B) is a preferred curvature and r\ a is a positive 
constant [30] • Qo, Qi and B must be obtained by fitting to experiments. Q\L 2 + 
2QoQiL cos(B) < means that the curvature amplitudes at the head are larger than those at 
the tail and we should expect traveling waves going from head to tail causing the nematode 
to swim forward; similarly, Q 2 L 2 + 2Q QiLcos(B) > should give traveling waves going 
from tail to head so that the nematodes swim backward \12\. The total moment M(s, t) can 
be written as 

M = M p + M a = EI(K - K a ) + (V P - Va)I~^ = EI(K - Ka) + 7? J — . (A.6) 

Note that if rj = 1] P — i] a > then there is net dissipation of energy in the tissue; if 
77 = T] p — T] a < then there is net generation of energy in the tissue. For live nematodes 
that actively swim in the fluid we expect 77 = rj p — i] a < since the net energy produced 
in the (muscle) tissue is needed to overcome the drag from the surrounding fluid. In other 
words, the driving force for the traveling waves seen in the nematode has its origins in the 
contractions of the muscle. 

We model the nematode as a hollow cylindrical shell with outer radius r Q and inner radius 
Ti [19] such that the principal moment of inertia (second moment of the area of cross-section) 
I along the entire length of the nematode is given by I = | [(r m + |) 4 — (r m — |) 4 ] where 
r m is the mean radius and t is the cuticle thickness. 

Consistent with experimental observations, we assume that the nematode moves along 
the x-axis and that the deflections of the centerline of the nematode from the x-axis are 
small. This allows us to take s = x and write n(x,t) = d<p/dx = d 2 y/dx 2 where y(x,t) is 
the deflection of the centerline of the nematode and 0(x, t) is the angle made by the tangent 
t to the x-axis. The equations of motion can then be written as 

°-ir = c ' v - ^ = c -% < A7 > 

ox ox at 

9M IK s 

-^ + F V = 0, (A.8) 
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where F = F x e x + F y e y . V x is the velocity of the nematode along the x-axis and corresponds 



to the average forward speed U. Combining a linearized formulation of Eqs. (A.l ) and (A. 2) 



along with the viscoelastic model of Eq. (A. 6) offers a direct route towards (i) a closed-form 



analytical solution for the curvature K,(x,t) and (ii) an estimate of tissue properties (E and 
77). Note that due to the assumption of small deflections [26J, the x-component of the force 
balance becomes decoupled from the rest of the equations. As a result, we will not be able 



to predict V x even if y(x,t) is determined. But, we can solve equations (A. 6), (A. 7), and 



(A. 8) for appropriate boundary conditions on F y and M to see if we get solutions that look 
like traveling waves whose amplitude is not a constant, but in fact, is decreasing from the 
head to the tail of the nematode (Fig. [2^). To do so, we observe that the nematode's body 
oscillates at a single frequency irrespective of position x along its centerline (Fig. |2}d). We 
therefore assume y(x,t) and M(x,t) to have a form that involves a single frequency w, so 
that y(x, t) = f(x) cosut + g(x) smut, and 

M(x,t) = K b (u)-^cos(ut + iP{uj)) 

+K b (u)-^sia{ut + ip{u)), (A.9) 

where f(x) and g(x) are as yet unknown functions, ip(uj) is frequency dependent phase 
angle, and K b {uj) is a frequency dependent bending modulus of the homogeneous viscoelastic 
material making up the nematode. In particular, 

tanV = ^, K b = I^/W+oShf, (A.10) 
E 

where the parameter K b , E, and 77 correspond to the effective tissue properties of the nema- 
tode. To be consistent with the observation of a single frequency u and non-zero amplitudes 
of the curvature k at the head and tail, we apply boundary conditions 

F y (0,t) = 0, M(0,t) = 0, F y (L,t)=0, M(L,t) = 0, (A. 11) 



We can make further progress by differentiating the balance of moments once with respect 
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to x and substituting the y-component of the balance of forces into it to get 

d 2 M dy 



Substituting for M(x,t) from eqn.(A.9) yields a biharmonic equation for the functions / 



and g which can be solved to give the following solution for y(x, t) 

TT ll) 71 lb 

y(x,t) = A 01 exp(Pxcos(- + -)) cos(/3xsin(- + — ) - ut - O i) 

o 4 o 4 

71 lb 7T lb 

+ A 23 exp(-/3x cos(- + -)) cos(-/5x sin(- + —) - ut - 2 3) 
o 4 o 4 

. ,ip 3rr . . ip 3n . 

+ A 45 exp(^xcos(- - — )) cos(/?a;sm(- - — ) - cot - </> 45 ) 
4 o 4 o 

lb 37T w 37T 

+ A 67 exp(-/?xcos(- - — ))cos(-/3zsin(- - — ) -ut-(b 67 ), (A.13) 

where /3 = (Cntu/Kb) 1 ^, and Aoi, A 2 3, A 45 , A 67 , O i, 023, 045 and 067 are eight constants to 
be determined from the eight equations resulting from the sin cut and cos ut coefficients of the 
boundary conditions. Note that these are four waves of the type y(x,t) = b(x) cos(2tt(x + 
V w t) I A) which is the form originally assumed by Gray and Hancock based on experiment 
[26] . In contrast, we have obtained such waves as a solution to the equations of motion. We 
plot the amplitude and phase of these waves for a particular choice of parameters in Fig. |6j 
But, the exact solution above is cumbersome to use. We need simple expressions that can 
be easily fit to some observable in the experiment to obtain (3 and ip, or equivalently, E and 
rj. One such parameter is the amplitude of the traveling waves as a function of position x 
along the nematode. We develop a strategy to obtain this amplitude in the following. 

Based on the exact solutions to the equations we approximate the displacement y(x, t) 
as follows: 

/ \ a in /ip 3tt . . .ib 3tt 

y(x,t) W A\ exp(— px cos(- — )) cos(— px sm(- — ) — out) 

4 8 4 8 

lb 7T lb 71 

+ A 2 exp(/5x cos(— + -)) cos(Px sin(— + — ) - ut + 2 ) 
4 8 4 8 

yj 3tt ip 3tt 3tt 

+ A 2 exp((3x cos(- - — )) cos((3x sin(- - —) - ut + (b 2 - —) 

+ Aiexp(-/3xcos(| + |))cos(-/3xsin(^ + |)-^-|). (A.14) 
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For experimental values of K a (L, t)/n a (0, t) = 0.33 we find A\jAi ~ 443 so that A2 is much 
smaller in comparison to A 1 . Note that 

3"7T 7T 

?/(0, t) = A\ cos + A 2 cos(ut — 2 ) + A 2 cos(u;t — 2 + — ) + A\ cos(c<jt + — ) 

= v^Aicos^t + -) + 2A 2 cos( — )cos(cjt - 02 + — ). (A. 15) 

4 8 8 

The curvature can be calculated as under: 

. . d 2 y 2 . ,ib 37r . -0 37r 3n %b 

K(x,t) = -— ~ p Ai expl — px cos! )) cosi— px sin ) — ut 1 — ) 

v ' ; dx 2 PV H V 4 8 JJ V H V 4 8 ; 4 2 J 

+ /3 2 A 2 exp(pxcos(- + -)) cos(pxsin(- + -) - + 2 + - + -) 

n2 a in /"0 37T . /0 37T 37T 37T "0 

+ /3 A 2 exp(pxcos(- - y))cos(pxsm(- - — ) - + 2 - — - — + -) 
+ /3 2 v4 ie xp(-/5xcos(^ + -)) C os(-#rsin(J + + - + |)(A.16) 

Note again that 

re(0, t) = p 2 A 1 cos(cut + y - ^) + /5 2 A cos(cut - 2 - | - |) 

+/5 2 A 2 cos(a;t -^2 + ^- + ^--^)+ /5 2 ^i cos(c;t + | - |) 

= v^/S 2 ^! cos(a;t + - - J) + 2A 2 cos(— ) cos(c;t - 2 - j- + — ). (A.17) 
2 2 8 2 8 



It is possible to determine (3 and ib from (A. 15) and (A.17) alone. If we recognize that 



A\ » A 2 at the head then we see that the ratio of the amplitude of the curvature to the 
amplitude of the displacement at the head is simply f3 2 and the phase difference between 
them is (f — f ) (the phase difference between ^ and at the head is ^ — j). We find 
by comparing with the exact solution that an estimate of j3 using this method is accurate 
to within 1% and that of ib is accurate to within 2 or 3 degrees. 

The curvature is an oscillatory function and we can determine its amplitude A(x) simply 
by isolating the coefficients of cos cut and sin cut and then squaring and adding them. We get 
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the following expression as a result of this exercise: 

A 2 (x) = p A A 2 1 exp(-2Pxcos(^- 3 ^))+p 4 A 2 2 exp(2pxcos(^ + ^)) 

+/3 4 ^exp(2/3a;cos(^ - ^)) + ^A\ exp(-2/3xcos(^ + £)) 

-2/5 4 A 1 A 2 exp(-v / 2 / 9xsin(^ - J)) cos(-v / 2 / 9x sin(^ - J) - 2 ) 

4 8 4 8 

+2/3 4 A 1 A 2 cos(-2/3xsin(^ - ^) - 2 + ^) 

4 8 4 

-2/^|exp(V2/fccos(| - |)) cos(V2/fc cos(^ - |) + ^) 

+2/3 4 A 1 A 2 cos(2/3x sin(^ + + J + 2 ) 

4 8 2 

+2/5 4 A 1 A 2 exp( v / 2 / 9xsin(^ - ^)) cos(v / 2 / 9x sin(^ - + 2 - ^) 

4 8 4 8 4 

+2/3 4 ^exp(- v / 2>rcos(^ - ^))cos(v^/3xcos(^ - -) - -). (A.18) 

4 8 4 8 2 

We can fit the experimental data of curvature as a function of x using the expression above. 
There are five fit parameters - A 1 , A 2 , (3, ip and 2 . The value of 2 mostly affects the 
curvature profile near the tail. Values of 2 ~ 7r/3 seem to give good fits for the curvature 
data of the nematodes. A\ can be determined from the amplitude of the displacement at 
the head. This leaves three fit parameters - A 2 , (3 and ip. We can use this fit to check if the 
parameters j3 and tp obtained from analyzing the motion of the head alone are reasonable 
or not. 



Methods Summary 

C. elegans strain. All strains were maintained using standard culture methods and fed 
with the E. coli strain OP50. The following muscular dystrophic (MD) strains were used: 
fer-l(hc24ts), dys-l(cxl8)I and hlh-l(cc561)II;dys-l(cxl8)I double mutant. Note that hlh-1 
is a myoD mutant that qualitatively reveals the motility defects of dys-1 mutants. Analysis 
were performed on hypochlorite synchronized young adult animals, fer-1 mutants were 
hatched at the restrictive temperature of 25°C and grown until they reach the young adult 
stage, dys-1 (cxl8) /; hlh-1 (cc561) II mutants were grown at the permissive temperature 
of 16°C. Wild-type nematodes grown at the appropriate temperature were used as controls. 
Strains were obtained from the Caenorhabditis elegans Genetic Stock Center. 
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FIGURE LEGENDS 

Figure [TJ 

Motility of wild-type C. elegans and dys-l;hlh-l muscular dystrophic (MD) 
mutant swimming at low Reynolds number, (a) and (c): Visualization of C. elegans 
motion illustrating instantaneous body centerline or skeleton. Also shown are the nematode's 
(i) centroid and (ii) tail-tip trajectories over multiple body bending cycles, (b) and (d): 
Color-coded temporal evolution of C. elegans skeletons over one beating cycle. Results 
reveal a well-defined envelope of elongated body shapes with a wavelength corresponding 
approximately to the nematode's body length. 

Figure [2} 

Spatio-temporal kinematics of C. elegans forward swimming gait, (a) Represen- 
tative contour plot of the experimentally measured curvature {k) along the nematode's body 
centerline for approximately 6 bending cycles. Red and blue colors represent positive and 
negative k values, respectively. The y-axis corresponds to the dimensionless position s/L 
along the C. elegans' body length where s = is the head and s = L is the tail, (b) 
Nematode's body bending frequency obtained from Fast Fourier Transform of k at different 
s/L. The peak is seen at a single frequency (~ 2.4 Hz) irrespective of the location s/L. 
(c) Contour plot of curvature k values obtained from the model. The model captures the 
longitudinal bending wave with decaying magnitude, which travels from head to tail, (d) 
Comparison between experimental and theoretical curves of k at s/L = 0.1 and s/L = 0.4; 
dashed lines correspond to model predictions (root mean square error is ~ 10% of peak-to- 
peak amplitude). 

Figure [3f 

Schematic of the analytical model for the total internal moment M. Muscle tissue 
is described by a visco-elastic model containing both passive and active elements. The 
passive moment (M p ) is described by the Voigt model consisting of a passive elastic element 
(spring) of stiffness E (i.e. Young's modulus) and a passive viscous element (dashpot) of 
tissue viscosity r] p . The active moment (M a ) is described by an active muscular element of 
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viscosity rj a and illustrates a negative slope on a force- velocity plot. Since there is a net 
generation of energy in the muscle to overcome drag from the surrounding fluid, we expect 
V = V P ~ Va < 0. 

Figure [4f 

Typical C. elegans viscoelastic material properties. Typical experimental profile of 
the curvature amplitude \k\ decay as a function of body position s/L. Color-coded theoretical 
profiles of the curvature amplitude decay \k\ at fixed value of the bending modulus K b . 
Curves vary from ip = 0° (red) to ij) — —90° (blue), which corresponds to i] = and E — 0, 
respectively. 

Figure [5} 

Kinematics and material properties of wild type and three muscle mutants of C. 
elegans. (a) Measured kinematic data and (b) estimated Young's modulus E and absolute 
values of tissue viscosity \rj\ for wild-type, fer-l(hc24), dys-l(cxl8), and dys-l(cxl8); hlh- 
l(cc561) adult nematodes (n = 7-25 nematodes for each genotype. * - p < 0.01). 

Figure [6j 

Amplitude and phase of traveling waves obtained from enforcing force and mo- 
ment boundary conditions at x = 0, L. (a) Amplitudes of the A 2 3 and A e7 waves 
decrease from x = to x = L while the amplitudes of the A i and waves increase from 
x = to x = L. The amplitudes of the former two waves are larger than the latter two. 
(b) Phases 0oi, 023, 045 and 067 are all constant from x = to x = L. The phase difference 
between the Aqj wave and A23 wave is approximately 90°. The parameters used to obtain 
these plots are: u = Air radians/sec, L = 1.0mm, Kf, = 5.0 x 10~ 16 Nm 2 , C n = 0.06Ns/m 2 , 
i) = -45°, EIQ = 4.35 x lQ- 12 Nm, Q X L = -1.054Q and B = 198.4°. 
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